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ABSTRACT 


The  paper  is  a  survey  of  some  aspects  of  information  based  selection  of 
numerical  methods.  Some  theoretical  ideas  on  how  to  deal  with  uncertain 
information  are  discussed  and  the  example  of  an  universal  quadrature  formula 
is  introduced.  Certain  aspects  of  the  selection  of  finite  element  method  and 
adaptive  mesh  construction  are  discussed.  Numerical  examples  illustrate  the 
theoretical  aspects.  ^  '  ' 
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INTRODUCTION 


Any  practical  computation  is  aimed  to  obtain  reliable  results  in  an 
optimal  way.  The  notion  of  the  optimality  is  very  complex  in  practice.  It 
includes  not  only  computer  cost  (for  example,  depending  on  various  aspect  of 
hardware),  but  includes  also  the  manpower  cost  (salaries).  Manpower  cost 
makes  today  typically  90-95?  of  the  total  cost  of  an  engineering  project.  The 
optimal  choice-optimization  of  the  method- strongly  depends  on  the  goal  of  the 
computation  and  on  the  available  information.  An  essential  part  of  the 
computation  is  also  the  assessment  of  the  reliability  of  the  obtained 
results.  Let  us  refer  to  [Noor,  Babuska  (1987)]  which  gives  a  survey  of  the 
principles  of  the  quality  assessment  of  the  finite  element  computation  and 
major  literature  (about  200  references  cited).  This  survey  addresses  many 
aspects  inclusive  adaptive  approaches,  optimal  extraction  of  the  desired 
information  from  the  computed  data  (in  the  post  processing  phase),  a- 
posteriori  error  analysis,  etc. 

In  practice  any  computation  is  information  based.  Many  numerical 
methods  and  codes  for  solving  mathematically  formulated  problems  are  available 
today.  A  very  important  problem  is  to  characterize  the  co-'ditions  under  which 
a  concrete  method  performs  well.  In  engineering  a  large  effort  is  spent  to 
obtain  indirectly  such  characterization  by  comparative  computational  studies 
of  various  benchmark  problems.  See  e.g.  [MacNeal,  Harder  (1985)],  [Robinson, 
Blackham  (1981)]  any  many  others.  Unfortunately,  mathematical  theory  here  is 
practically  nonexistent. 

Presently  the  research  is  also  going  on  how  to  use  the  expert  systems 
(in  general,  application  of  the  principles  of  artificial  intelligence)  to 
select  a  numerical  method,  its  part  or  other  basic  parameters  of  the 

V 

computational  analysis.  See  e.g.  [Babuska,  Rank  (1987)],  [Rogers,  Barthelemy 
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(1985)],  [Dyra  (198M)]  and  others. 

Mathematically  the  optimal  selection  of  numerical  methods  was  for  long 
time  under  consideration.  See  e.g.  [Babuska,  Sobolev  (1965)],  [Bachvalov 
(1968)],  [Sobolev  (197^)],  [Traub,  Wozniakowski  (1985)],  for  basic  ideas, 
results  and  literature.  We  will  address  some  aspects  of  the  optimization  of 
numerical  methods  below. 

Information  about  the  class  of  problems  (or  set  of  their  solutions) 
under  consideration  is  essential  for  the  optimal  selection  of  a  method.  For 
theoretical  aspects  of  the  notion  of  an  information, optimal ity  and  information 
based  complexity,  we  refer  to  [Wozniakowski  (1985)]. 

It  is  imperative  to  concentrate  on  information  which  is  practically 
available  and  realistic  and  not  only  mathematically  elegant  and  convenient. 

The  tests  on  benchmark  model  problem  reflecting  the  practice  are  necessary  to 
keep  the  research  in  the  prospective.  The  optimal  method  selection  depends 
strongly  on  the  set  of  solutions  (information).  Some  mathematically  elegant 
results  could  be  practically  misleading  if  taken  out  of  context.  For  example 
[Smoljak  (1966)]  has  shown,  in  the  case  of  quadrature  formulae,  that  for 
convex  balanced  set  of  Integrated  fiinctions  the  linear  algorithm  is  as  good  as 
nonlinear  (adaptive)  one.  Nevertheless,  he  (see  also  Bachvalov  [1968]) 
underlined  that  result  of  this  type  is  very  information  dependent  and  has  not 
to  be  overstated  in  general. 

Heuristics  is,  and  has  to  be,  used  directly  or  indirectly  in  the 
selection  of  an  optimal  method.  This  heuristics  is  either  in  the  selection  of 
the  available  information  (although  following  mathematical  theory  can  be 
rigorous)  or  in  the  optimization  reasoning  itself  (and  it  Is  not 
mathematically  rigorous).  A  very  elegant  analysis  based  on  a  probabilistic 
justification  of  practical  heuristic  arguments  is  in  [Gao  (1986)]. 


'  ^  •  «  F  *  F  •  •  m  *  *  •  •  •  •  '  »  •  •  ■  m  • 
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We  will  address  here  some  problems  and  mathematical  results  and  will 
present  some  basic  ideas  on  sample  examples.  For  simplicity  we  will  discuss 
some  aspects  related  to  3  directions. 

a)  Selection  of  nonadaptive  methods. 

b)  Feedback  and  adaptive  methods. 

c)  Reliability  assessment. 
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2.  SELECTION  PRINCIPLES  FOR  NONADAPTIVE  METHODS 
2.1.  THE  QUADRATURE  FORMULA 

The  simplest  example  studied  in  detail  in  the  literature  is  the  problem 
of  the  optimal  quadrature  formula.  This  example  can  also  serve  as  a  prototype 
of  various  approaches  used  in  more  complex  setting. 

The  notion  of  the  optimality  is  very  broad  and  relates  to  different 
types  of  convergence.  See,  e.g.  [Sobolev  (1962a)],  [Sobolev  (1962b)]  and 
others.  For  an  extensive  theory  of  optimal  cubature  formulae  and  the 
functional  analytic  prerequisites  for  the  study  of  optimal  formulae,  we  refer 
to  the  large  monography  (808  p)  of  [Sobolev  (197^1)]. 

Let  us  address  some  aspects  of  quadrative  formula  in  its  most  elementary 
setting.  Let 

F(u)  -  J  u(t)dt 
0 

be  a  functional  defined  on  a  certain  space  of  continuous  functions  on  I  » 
(0,2ir).  Assume  that  the  values  u(tj)  are  computable  for  any  0  <  tj^  <  2it 
(and  only  these  values  are  computable).  Let  us  be  Interested  in  the 
quadrature  formula  using  N  function  values  with  the  minimal  error.  This 
leads  to  the  following  problem.  Given  N  >  0,  integer,  find  0  ^  t^  i  2it, 


,  i  »  1,2,...,N  such  that 

N 

(1)  |F(u)  -  I  u(tj^)a.|  »  E(u,N) 

i»1 


is  minimal  over  a  set  S  of  functions  u,  i.e.,  let  us  be  interested  in 

(2)  inf  sup  E (u,N) . 

t^,a^  u€S 


The  formula  which  achieves  (2)  will  be  called  optimal  and  denoted  by  Q(S,n). 


Essential  here  is  the  set  S.  It  is  typically  selected  as  the  unit  ball  in  a 
Banach  space.  In  (2)  we  have  taken  the  infimum  over  all  0  ^  tj^  ^  2ii  and 
without  any  constraints.  We  can  be  also  interested  in  the  case  when  t^ 
are  constrained  to  uniform  mesh,  etc.  For  various  results  related  to  the 
optimality  of  this  type,  we  refer  to  [Sobolev  (197^)],  Nikolski  (1958)], 
[Traub,  Wozniakowski  (1980)]  and  others.  The  main  problem  for  the  application 
of  the  selection  of  a  formula  based  on  this  (worst  case)  optimality  principle 
is  the  selection  of  the  set  5.  To  illustrate  this  difficulty,  let  us 
consider  one  parametric  family  (scale)  of  spaces 

5k  *  I  M  k  - 
h‘^(I) 

where  h'^(I),  k  >  ‘4  is  the  standard  Sobolev  space.  The  optimal  formula 
Q(S|^,N)  depends  on  k  and  N.  Hence,  the  following  natural  question  is: 

For  a  concrete  u,  what  k  to  select?  In  this  connection  we  can  deal  also 
with  countably  normed  spaces,  i.e.  to  assume  that 

|u|  i  iKk),  *4  <  k  <  », 

with  a-priori  given  "Kk).  For  more  we  refer  once  more  to  [Sobolev  (197^)] 
and  also  to  [Babuska,  Prager,  Vitasek  (1966)]. 

It  seems  to  be,  that  the  selection  of  the  quadrature  formula  as  an 
optimal  one  with  respect  to  a  special  choice  S|^,  is  practically  ineffective 
because  of  the  uncertainty  In  the  selection  of  Sj^. 

Let  us  discuss  this  aspect  on  the  problem  when  the  function  u  is  2Tr- 
periodlc.  The  problem  of  integration  of  a  periodic  function  was  analyzed  in 
many  papers.  For  example  [Bachvalov  (196il)]  studied  optimal  lower  and  upper 

(r^  ,  .  •  .  ifg) 

bounds  for  the  errors  in  s-dimensional  cases  of  classes  H  of 

s,p 


functions  u(x^,...,Xg)  having  period  p  in  every  direction  and 

(I'l  »  •  •  • 

|f  II  <  1 •  The  estimates  were  the  best  possible  ones  up  to  the 

term  log  N.  The  upper  estimates  were  obtained  by  number  theoretical 
approaches  of  [Korobov  (1963)]  and  the  lower  ones  were  studied  by  the  theory 
presented  in  [Kolmogorov,  Tichomirov  [1959)].  Many  other  results  are 
available. 

Let  us  return  now  to  the  one  dimensional  case  of  periodic  function  and 
address  the  problem  of  uncertainty  of  the  apace  selection.  For  a  detailed 

V 

theory,  see  [Babuaka  (1968)].  First  we  could  ask:  What  is  the  intuitive 
content  of  the  statement  that  a  (complex)  function  u  is  a  2ir-periodic , 
continuous  one?  We  can  formulate  it  so, that  it  belongs  to  a  Hilbert  space 
H  which  has  the  following  property  P: 

:  If  f  €  H,  then  f  is  continuous. 

P2:  H  is  dense  in  C. 

:  If  c  is  real,  f  €  H,  then  also  g(x)  =»  f(x'*-c)  €  H  and 

P^ :  There  exists  K(H)  such  that 

|f|c  <  K(H)|f|j^  . 

(We  restrict  ourself  on  the  case  of  Hilbert  spaces  only). 

We  have 

THEOREM  1 .  Let  H  be  a  Hilbert  space  having  property  P.  Then 
(3a)  1)  €  H,  k  -  ...,-1,0,1,... 

(3b)  2)  functions  k  -  ...-1,0,1...  create  orthogonal  basis  of  H 

(3c)  3)  denoting  n|^  =•  .  then 


Space  satisfying  the  condition  P,  respectively  (3a),  (3b)  and  (3c)  will  be 
called  a  periodic  space.  We  will  need  also  a  stronger  notion.  The  space  H 


is  called  strongly 

periodic, 

if  it  is  periodic  and 

(3d) 

■  il<X| 

1  “ikX| 

1«  l„  ■ 

|e  Ih 

(3e) 

If  |j|  i 

|k|,  then 

4-00 

(3f) 

t=0 

^  0  ^  a  ^  2 

n  >  1  integer 


(6) 

where 

(7) 


A(n,H)  =  sup  |F(u)  -  Tj^(u)| 

M„<' 

J’’ 


is  the  trapezoid  formula.  Obviously  u  is  related  to  the  optimal  formula  for 
all  distribiitions  of  tj ,  6  when  we  restrict  ourself  to  uniform  mesh  and 
A  is  the  error  for  the  trapezoid  formula. 

Intuitively,  one  can  expect  that  the  trapezoid  formula  compares  well 


with  the  optimal  one.  Nevertheless,  we  have 


arbitrary.  Then  there  exists  a 


THEOREM  2.  Let  >  0,  J  -  1,2... 
periodic  space  H  such  that 


(8) 


lim  3UP 


THEOREM  3.  Let  H  be  any  strongly  periodic  space.  Then 


(9) 


A(n.H) 

lim  sup  - - -  < 

n-*®  aj(n,H)/n 


THEOREM  4.  For  every  period  space 
(10)  A(n,H)  >  6(n,H). 

THEOREM  5.  Let  H  be  a  periodic  space.  Then 


(11) 


lim 

n-+® 


A(n.H) 

6(n,H) 


1  . 


Only  the  trapezoid  formula  has  the  property  that 


(12) 

for  all  periodic  spaces. 


lim  sup 
n-*^“ 


A(n.H) 

5(n,H) 


< 


CD 


These  theorems  show 

1)  Speaking  about  periodic  function  u,  we  likely  mean  that 
to  a  strongly  periodic  space. 

2)  Although  the  trapezoid  formula  is  not  optimal  ‘"or  any  H, 
good  for  all  periodic  spaces,  i.e.  only  the  trapezoil  fo’^nula  is  a 
f  ormula . 

The  theorems  show  that  the  selection  of  an  optimal 
particular  space  is  likely  a  bad  choice  in  practice. 


□ 

□ 


n 

it  belongs 

it  is 
robust. 


*■  ern  ila  for  a 
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Let  us  illustrate  it  in  a  concrete  example.  Let  |u|„  - 

2it  ” 

J  (u^  +  u'^)dx,  then  the  formula 
0 

(13)  C(n,H)T^(u) 

(14)  C"'’(n,H) 

is  the  optimal  one,  i.e. 

sup  |F(u)  -  C(n,H)T  (u)|  »  p(n,H). 

|u|„0 

Assume  now  that  u  =  e“  a  -  3,10.  Then 

^  F(e^  -  4.88079258586502408 

iiir 

1  p^gio  sin  x^  ^  2815.71662846625447 

cir 

and  the  values  obtained  by  the  trapezoid  and  the  optimal  formula  are  given  in 
Table  1 . 


TABLE  1 


TRAPEZOID 

FORMULA 

OPTIMAL  FORMULA  j 

n 

a  -  3 

a  »  10 

a  •  3 

a  -  1  0 

8 

4.88241999058958100 

2047.90959481  962441 

4.64604604. . 

2900.35030. . 

16 

4.880792585936661  73 

2815.77672896656761 

4.81 902223.. 

2780.14081  . ., 

'24 

{ 

4.88079258586502408 

2815.71662897903758 

4.85310536. . 

2799.74394. . 

’  J 


1 


n^  t.1  (tA)^  +  (1^) 
n 


Table  1  suggests  that  in  practice  only  for  small  n  and  small  accuracy 
the  "optimal"  formula  can  be  better,  provided  that  the  space  H  is  properly 
selected  (which  will  not  likely  to  happen).  For  detailed  theory  and  many 
additional  aspects,  we  refer  to  [Babuska  (1968)]. 


The  observation  about  the  universal  property  of  the  trapezoid  formula 
and  its  analysis  was  made  later  in  various  settings  by  various  authors.  See, 
e.g.,  [MotornyJ  (197^)]  and  others.  We  refer  here  also  to  [Traub, Wozniakowski 
(1980)]  where  many  papers  about  optimal  formulae  are  cited. 

Similar  principles  of  (worst  case)  optimality  can  be  applied  in  many 
more  complex  cases.  We  mentioned  the  case  of  optimal  difference  formula,  the 
optimal  choice  of  the  trial  functions  in  variational  methods,  etc.,  see,  e.g. 
[Babuska,  Prager,  Vltasek  (1966)],  solution  of  initial  value  problem  for 
ordinary  differential  equation,  see,  e.g.  [Bachvalov  (1963)],  integral 
equations,  see  [Emeljanov,  II Jin  (1967)]  and  many  others. 

2.2.  THE  REGULARITY  OF  THE  SOLUTION  AS  THE  BASIC  INFORMATION  FOR  FINITE 
ELEMENT  SELECTION 
Let  us  consider  the  model  problem 


(15a) 

-Au  • 

f 

on 

a. 

(15b) 

u  » 

g 

on 

da. 

where  c  R^ 

is  a  bounded  domain 

and 

da 

is  its  boundary. 

The  finite  element  method  consists  (in  the  most  simple  case)  in  the  par¬ 
tition  of  q  into  the  set  of  triangle  and  quadrilateral  elements  (the 

mesh)  and  in  the  best  approximation  of  u  (in  H^(q))  by  piecewise  poly¬ 
nomials  of  degree  p^^  on  t^.  The  space  of  these  piecewise  polynomials  is 
called  trial  space  and  its  dimension  N  is  called  the  number  of  degrees  of 
freedom. 

Typically  it  is  assiimed  that  the  only  Information  available  about  u  is 
that  it  belongs  to  H^(q),  k  >  1  and  the  optimal  error  (for  a  quasiuniform 


mesh)  is 
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(I6a) 


Nl  1  <  CN  |u| 

H  (£5)  H'^(n) 


where 


(16. b) 


li  “  mln(p,k-1) 


when  P|  »  p  for  all  elements  This  error  Is  the  optimal  one.  Removing 

k-1 

constraint  Pj^  =  Pt  the  best  estimate  is  still  o(n  )  which  can  be 
easily  be  proven  by  concept  of  the  n-width. 

If  no  other  information  than  u  €  H*^(£J)  is  available,  then  the  uniform 
partition  is  obviously  preferable  because  it  leads  to  the  best  possible 
estimate.  For  more  details,  see  [Babuska,  Suri  (1987)]. 

The  Information  that  u  €  H*^(P)  is  very  far  from  the  optimal  one  in 
practice  (e.g.  structural  mechanics).  Usually  the  data  (l.e.  boundary,  f) 
are  piecewise  analytic.  Then  the  solution  of  the  problem  (15a, b)  belongs  to  a 
countably  normed  space.  It  is  possible  to  prove  that 


(17a) 


J  I  D*^u|^<Ji^_2,.g(x)dx  i  Ck!d*^, 


k  -  2,  . , 


where 


(17b) 


|k-2*B 


and  are  typically  the  vertices  of  £1. 

For  the  proof,  see  [Babuska,  Guo  (1985)],  [Babuska,  Guo  (1986)].  We 
have  now 


THEOREM  6.  Let  u  f  hV  £1)  satisfies  (’7).  Then  there  e'^ists  sequence 


of  meshes  and  elements  degrees  on  such  that 
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(18)  |e|  ,  <  Ce"^’^ 

h'(Q) 

where  Y  >  0  depends  on  0  and  d  in  (17).  □ 

The  exponential  rate  (18)  can  be  obtained  in  practical  computations. 

For  more,  see  [Guo,  Babuska  (1986a)],  [BabuSka,  Guo  [1986)].  Results  of  this 
type  were  implemented  in  the  commercial  finite  element  code  PROBE  (Noetic 
Tech.,  St  Louis),  [Szabo  (1985)].  For  the  survey  paper  about  the  state  of  the 
art  of  the  h-p  version  of  the  finite  element  method,  we  refer  to  [Babuska 
(1986)]. 

The  decision  in  which  space  the  solution  u  should  be  imbedded  is 
crucial  for  selection  of  the  finite  element  meshes  and  degrees.  In  Fig.  1  we 
show  the  accuracy  of  the  computation  of  the  elasticity  problem  of  a  L-shaped 
domain  by  various  finite  element  methods  (meshes,  degrees)  as  function  of  N 
(which  roughly  expresses  the  cost)  and  which  are  optimal  for  various  selection 
of  the  spaces  of  the  solutions.  This  shows  similar  aspects  we  have  addressed 
in  previous  section,  namely  very  different  performances  of  the  method. 

In  one  dimensional  setting  much  more  details  are  available.  See  [Gui, 
Babuska  (1986)].  Consider  the  following  sample  problem:  Let  u(x)  »  x“  - 
X,  a  >  ‘4  be  given  on  I  -  (0,1).  Consider  the  set  of  all  partitions  of  A 
of  I 

A  0  >  Aq  <  ...  <  A„  -  1,  Ij'  - 

Denote  by  S(A,£)  -  (u  €  Hq(I)|u  is  polynomial  of  degree  p^  on  I^}. 

Let  us  be  interested  now  on  what  can  be  said  about 
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dim  S( A,£)*N 


We  have  then 


THEOREM  7. 


c,(„,  _L_g^(^  <  Q  < 

/iF^ 


2 

where  gQ  -  (/2  -  1 )  . 
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Fig.  1  .  The  error  In  the  energy  norm  In  dependence  on  N  for  various  meshes 
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THE  ADAPTIVE  HETHODS 


It  is  worthwhile  to  distinguish  between  a  feedback  method  and  an 
adaptive  method.  A  feedback  method  is  any  method  which  utilizes  the  computed 
values  to  steer  its  direction.  An  adaptive  method  is  a  feedback  method  which 
is  optimal  in  a  precise  sense.  For  more  elaboration,  we  refer  to  [Rheinboldt 
(1983)].  [Babuska  (1986b)]. 

3.1.  THE  FINITE  ELEMENT  METHOD 

Let  us  explain  the  main  ideas  on  the  example  of  an  adaptive  finite 
element  solver  for  one  dimensional  boundary  value  problem 

(20a)  -(au’)  '  ♦  bu  -  f  x  e  I  -  (0,1 ) 

(20b)  u(0)  -  u(1)  -  0. 

The  finite  element  method  consists  of  the  mesh 

A  :>  O-Xq  <  X,  <  •••  <  Xf,  •  1 

of  the  cardinality  N  -  ic(A).  Denote  by  Ij  -  (xj_^,Xj)  the  element.  In 

the  most  simple  case  the  finite  element  solution  u(A)  c  Hq(I)  is  a 
continuous  function  on  I,  linear  on  Ij,  J  »  1 , . . . , N.  Given  the  mesh  A, 

the  finite  element  solution  u(A)  Is  uniquely  defined.  Denoting  by  Uq  the 

(exact)  solution  of  (20),  the  error  is  e(A)  -  u(A)  -  Uq.  Assume  that  the 
error  is  measured  In  the  energy  norm  l®|£  * 

|e|p  -  J  (ae'^  ♦  be^)dx. 

0 

The  feedback  method  now  consists  of  the  construction  of  a  sequence 
Ajt  J  •  1,2,...,  so  that  Aj  depends  on  Aj  and  u(A|^).  2.  <  J .  The 

operator  'Y  giving  Aj  is  called  transition  operator  and  defines  the 


feedback.  Sequence  {A^}  is  called  the  trajectory.  A  feedback  method  is 
called  adaptive  if  it  creates  an  optimal  trajectory.  Various  optimaliiy 
definitions  can  be  considered.  See  [Babuska  (1986b)] 

We  will  call  the  trajectory  {A^}  optimal  with  respect  to  the 
convergence  if 

(21  )  |e(Aj)  Ig  ■*  «  as  J  -►  ». 

Define  now 

(22)  ®,(N)  •  inf  |e(A)L 

•c(A)-N 

i.e.  ♦u(N)  is  the  smallest  error  which  can  be  obtained  by  the  mesh  with  the 
cardinality  N, 

We  will  call  the  trajectory  {A^}  optimal  with  respect  to  convergence 
rate  if 


(23) 


lim  sup 

i-»<» 


4  («(4  )) 


<  ». 


For  a  transition  operator  ct  defining  the  feedback  method,  we  define 
S(,(uq,A^  , a)  (respectively  Sr(uq, A^  ,rt))  -  {(uQ,A^,tt)|  the  trajectory  is  optimal 
with  respect  to  the  convergence  (convergence  rate)}.  The  goal  is  to 
design  ec  so  that  and  Sp  will  be  so  large  that  they  Include  all  cases 

which  are  important  in  practice.  In  [Babuska,  Vogelius  (1984)]  we  considered 
some  feedback  methods  (the  details  are  technical  and  cannot  be  given  here) 
which  have  the  following  properties: 

1)  Sp(uQ,  Ai  ,i''J)  =•  hJ(I).  that  the  feedback  leads  always  to  the 

convergence. 

2)  Sp(uQ,A^,tt)  y  Hg(I).  Sufficient  condition  for  Uq  ^  Sp  are 


given.  These  conditions  are  satisfied  for  solutions  which  are  important  in 


application  as,  e.g.,  Uq  »  x“  -  x,  a  >  *4  and  many  other  functions. 

Functions  u€Hq(I),  u  i  SR(uQ,A^,£t)  were  also  constructed.  For 

mathematical  details  see  [Babuska,  Vogelius  (198M)]. 

The  selection  of  the  method  depends  very  much  on  the  information  about 

the  solution.  This  information  is  usually  not  too  reliable  and  hence  robust 

methods  have  to  be  preferred;  in  our  case,  robustness  is  directly  related  to 

the  size  of  Sj^.  The  decision  about  the  available  information  will  always  be 

left  to  the  user.  Likely,  the  physical-engineering  arguments  are  the  most 

reliable  way  to  characterize  available  information.  Probability  approach  does 

not  avoid  this  problem  because  the  information  about  the  probability  field  is 

not  available  and  is  usually  made  on  very  arbitrary  basis. 

We  mentioned  one  dimensional  problem.  The  ideas  were  extended  to  more 

dimensional  problems.  See,  e.g.  [Babuska,  Miller  (1987)].  An  adaptive  code 

FEARS  (for  elements  of  degree  p  *  1)  has  been  written  and  experimented  with. 

Although  in  more  dimensions  such  a  detailed  and  complete  analysis,  as  in  one 

dimension,  is  not  available,  the  numerical  experiments  indicate  that  the  coded 

method  Is  adaptive  with  respect  to  the  convergence  rate  for  a  very  broad  set 

of  solution  of  engineering  interest.  Fig.  2  shows  the  finite  element  mesh  for 

the  elasticity  problem  of  a  cracked  panel.  The  exact  solution  has  singularity 

at  the  tip  of  the  crack  and  uniform  mesh  would  give  the  rate  0(N  ")  while 

-*/ 

the  adaptive  solver  leads  to  the  rate  0(N  *)  (which  is  the  same  as  for  the 

-V 

smooth  solution).  As  coul  ^  be  seen  from  Table  2,  the  rate  0(N  *)  was 
practically  achieved.  Table  2  shows  also  the  effectivlty  index  0  of  the 
estimator  as  defined  in  the  next  Section. 


Figure  2.  The  meah  constructed  by  FEARS  for  N  »  617. 


TABLE  2 


MESH 

N 

I'U  , 
We 

0 

1 

101 

26.38 

.885 

2 

143 

21 .35 

.991 

3 

221 

16.79 

1  .058 

i4 

301 

1  3.61 

1.116 

5 

61  7 

9.63 

1  .088 

As  in  one  dimension,  the  set  Sp  ^  H^(n)  and  hence  there  exists  cases  where 
the  method  does  not  lead  to  an  optimal  (convergence  rate)  trajectory. 
Nevertheless,  these  cases  are  very  likely  without  any  practical  importance. 
For  the  principles  of  adaptive  method,  we  refer  also  to  [Babuska,  Gui  (1986)] 
and  [Noor,  Babuska  (1987)]. 
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3.2.  THE  ADAPTIVE  ODE ’3  SOLVERS 

The  modern  standard  ODE's  solvers  are  of  feedback  type.  They  are 
usually  based  on  "per  unit  step"  or  "per  step"  tolerance  criterion.  The 
tolerance  t  (which  is  an  input)  has  in  principle  two  purposes:  the  accuracy 
control  or  the  feedback  control.  We  will  discuss  here  briefly  the  feedback 
aspects.  Let 

(23)  X  »  f(t,x),  x(0)  =  Xq,  0  <  t  <  T 

be  the  problem  under  consideration.  Then  approximate  solution  5(t) 
satisfies 

|(t)  =»  f(t,C)  ♦  n(t),  t(0)  -  Xq,  0  <  t  <  T 

and  we  can  judge  the  quality  of  5(t)  by  ri(t)  and  define  dist(x(  t ) ,  5(  t ) ) 

*  l^lLp((0,T))-  Now  we  have 

THEOREM  8.  The  "per  unit  step"  (respectively  "per  step")  approach  is 
adaptive  with  respect  to  the  rate  of  convergence  for  |ri|j^  ((0  T))  ’ 
respectively  |n|^^  ^  measure,  for  a  large  set  of  f. 

In  practice  It  is  not  necessary  to  find  exactly  the  optimal  mesh  because 
the  effectivity  of  the  solution  is  not  too  sensitive  to  the  perturbances  of 
optimal  mesh. 


ki 
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A  POSTERIORI  ASSESSMENT  OF  THE  ACCURACY 


The  a-posteriori  assessment  of  the  accuracy  of  the  computation  is 
essential  but  very  delicate.  For  various  approaches  we  refer  once  more  to 
[Noor,  Babuska  (1987)]  and  literature  cited  there.  In  general,  for  the  a- 
posteriori  error  estimate  we  desire  that  the  effectivity  index  0  of  the 
estimator  E 

estimate 

0  ”  r - 

true  error 

has  the  following  properties: 

<  0  <  K2 

where  K-j  and  K2  are  independent  of  any  detailed  characteristics  of  the 
mesh,  solution,  etc. 

b)  0  -►  1  as  E  -►  0 

and  for  reasonable  accuracies  .9  <  0  <  1,1  (say),  for  reasonable  problems 
and  meshes. 

It  is  practically  important  that  0  is  close  to  one  for  practical 
reasons.  Table  2  has  shown  this  effectivity  index  for  the  cracked  panel 
problem.  As  before,  we  wish  to  characterize  the  set  of  the  solutions  and  the 
meshes  for  which  the  above  properties  hold  (i.e.  we  characterize  analogs  to 


2 


5.  CONCLUSIONS 

Any  numerical  analysis  is  information  based.  Hence,  mathematical 
analysis  is  very  important, provided  that  it  addresses  the  pertinent  problems 
of  the  computational  analysis.  Nevertheless,  various  mathematical  results 
could  be  practically  misleading  if  their  assumptions  are  not  confronted  with 
the  circumstances  of  practical  computations.  Hence,  it  is  imperative  to  apply 
the  obtained  theoretical  results  in  the  environment  of  practical  computations, 
made  comparative  numerical  studies  and  test  the  applicability  of  the 
theoretical  conclusion.  Otherwise  the  results  of  information  based  complexity 
would  have  no  impact  on  computational  analysis.  In  the  paper  we  did  try  to 
address  some  theoretical  aspects  in  the  relation  with  computational  practice 
in  a  most  simple  setting. 
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